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In a recent paper, Vakarchuk [T] approached the Dirac equation coupled 
to a vector Coulomb potential for a particle with an effective mass dependent 
on the position. He assumed that the effective mass has a form of a multipole 
expansion and took only the first two lowest terms into account. In effect, 
he considered a fermion in the background of a mixed vector-scalar Coulomb 
potential. 

Back in 1973 Soff [2] et al. found the analytic solution of the Dirac equa- 
tion with an arbitrary mixing of vector and scalar potentials. The solution 
for the spin and angular variables was expressed in terms of spinor spherical 
harmonies, also called spherical spinors, resulting from the coupling among 
two- component spinors and spherical harmònic functions. The radial equa- 
tions for the upper and lower components of the Dirac spinor were treated by 
the brute force power series expansion method. In a recent time the solution 
of this problem was used to speculate about the breaking of pseudospin sym- 
metry in heavy nuclei [3]. It is also worthwhile to mention that a pedagogical 
approach to the Dirac equation coupled to a mixed vector-scalar Coulomb 
potential is already crystallized in a textbook 0]. 

It seems that Vakarchuk did not know about the solution of the Dirac 
equation with a mixed vector-scalar Coulomb potential and, as a matter of 
fact, he presented a more elegant method to find out the analytic solution. 
Although it may seem strange, the author tried to map the Dirac equation 
into an effective Schròdinger equation for all the components of the Dirac 
spinor. Nevertheless, I am afraid that something might be by far wrong. 

In order to clarify my criticism, let me begin writing the Dirac equation 
with the effective mass and the Coulomb interaction as given in PQ: 




where a and e are constants. With ijj defined as 
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the author of Ref. [T] obtains 
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where n = r/r, a stands for the Pauli matrices and the 4x4 matrices (3' and 
f3" are defined as 





Instead of a, a 4x4 matrix whose block diagonal elements are Pauli matrices, 
a, and whose off-diagonal block elements are zero, should used. Furthermore, 
contrary to the statement of author of Ref. [Tj, Eq. (jïïj) does not have the 
form of a Schròdinger equation. This criticism is endorsed by observing that 
the "centrifugal barrier" term contains off-diagonal matrix elements which 
mix the upper and lower components of the quadrispinor ip. 
By introducing the operator 

À = - (a.t+h) + Z -a.n (mc 2 a(3" + e 2 (3') (5) 

where L is the angular momentum operator, Eq. (jïïj) in the spherical coordi- 
nate system is supposed to reduce to (Tj: 
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where R is the radial part of ip, h 2 l* (Z*+l) is the eigenvalue of the operator 
À (k+h) for some eigenfunction, of course, and 
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with the upper sign for j — l + 1/2 and the lower sign for j = l — 1/2. Note 
that A operates on the spin and angular variables. It is a generalization 
of the spin-orbit coupling operator K, part of the toolkit for studying the 
Dirac equation with spherically symmetric potentials (see, e.g., and 
has the spinor spherical harmònic as eigenfunction. The eigenfunction of the 
operator A is not specified in Ref. jTj. In that case, what is the generalized 
spinor spherical harmònic? 

The author of Ref. [T| affirms that Eq. (jüj) formally coincides with the 
nonrelativistic Schròdinger equation for the Kepler problem. This decla- 
ration would be faithful if l* could only assume nonnegative integer vàl- 
ues but it does not in general. Therefore, one can neither identify R with 
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the usual radial functions of the nonrelativistic hydrogen problem R n i, with 
n = 0, 1, 2, . . . and l = 0, 1, 2, . . . , (n - 1) nor identify (E 2 - m 2 c 4 ) /(2mc 2 ) 
with the "Bohr formula" for the energy levels. Notice, though, that the 
identificat ion is believed to be true when e 2 = mc 2 |a|, i.e., when the scalar 
potential, if it is either attractive or repulsive, is as strong as the vector 
potential. 

If one insists that (jïïj) is the nonrelativistic Schrodinger equation for the 
Kepler problem then the Dirac eigenenergy solution (Eq. (4.2) in Ref. PQ) 
should have two branches of solutions, corresponding to positive and negative 
energy in a general circumstance. Furthermore, the condition to the existence 
of bounded solutions should be written as a < e 2 E/ (m 2 c 4 ) and not a < 
e 2 / (mc 2 ). 

Based on the above considerations it is not difficult to be vehemently 
opposed not only the solutions found in Ref. P but also the conclusions 
manifested there. 
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